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We present the first three-dimensional simulations of the evolution of a microquasar jet 
inside the binary-star system. The aim is to study the interaction of these jets with the 
stellar wind from a massive companion and the possible locations of high-energy emission 
sites. We have simulated two jets with different injection power in order to give a hint 
on the minimum power required for the jet to escape the system and become visible 
in larger scales. In the setup, we include a massive star wind filling the grid through 
which the jet evolves. We show that jets should have powers of the order of 10 37 erg/s or 
more in order not to be destroyed by the stellar wind. The jet-wind interaction results 
in regions in which high energy emission could be produced. These results imply the 
possible existence of a population of X-ray binaries not detected in the radio band due 
to jet disruption inside the region dominated by the stellar wind. 
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1. Introduction 

The powerful jets of X-ray binaries (XRB; i.e. microquasars) are produced close 
to a compact object, a black hole or a neutron star, through ejection of material 
accreted from the companion star. The interaction of jets with the stellar wind in 
the binary region has to be taken into account as a possible source of strong shocks 
and/or jet disruption in the case of a massive companion star. The occurrence of 
collisionless shocks can lead to efficient particle acceleration (see Ref. [I]), resulting 
in significant non-thermal emission of synchrotron and inverse Compton origin and, 
possibly, from proton-proton collisions (see Ref. 2 and references therein). 

In Ref. [3| (PBR08 from now on) the authors performed two-dimensional (2D) 
numerical simulations in order to show how the strong wind of an OB star can 
influence the jet dynamics at scales similar to the orbital separation (~ 0.2 AU). 
Simulations in two dimensions of a hydrodynamical jet interacting with an homo- 
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geneous (i.e. not clumpy) stellar wind were performed in cylindrical (axisymetric) 
and planar (slab) symmetry. The results showed that a strong recollimation shock 
is likely to occur at jet heights ~ 10 12 cm, which could lead to efficient particle 
acceleration and explain the TeV emission seen in several high-mass X-ray binaries 
(LS 5039, LS I +61 303, Cygnus X-l, see Refs. SJ [5]and[6]). It was also found that 
jet disruption could occur for jet kinetic luminosities as high as L } ~ 10 36 erg s _1 
because of jet instabilities produced by the act of the strong and assymetric wind. 
Such a Lj -value is ~ 0.1 — 1% the Eddington luminosity, typical for an X-ray bi- 
nary with persistent jets (see Fig. 1 in Ref. 7). This implies that the stellar wind 
in high mass microquasars could play a role not only for feeding accretion, but 
also for high-energy radiation and jet suppresion at the binary system scales. To 
carry this research further, we have performed three-dimensional (3D) simulations 
of hydrodynamical jets interacting with a strong stellar wind. We find, confirming 
previous results, that even jets with Lj > 10 36 erg s _1 may be disrupted, since the 
lower wind-jet momentum transfer in 3D as compared to 2D is balanced by the 
development of helical Kelvin-Helmholtz (KH) instabilities. Our findings also point 
to the presence of a strong recollimation shock that could efficiently accelerate very 
energetic particles. 

2. Numerical Simulations 
2.1. Set-up 

We have performed 3D simulations of two fast supersonic hydrodynamical jets with 
Lj = 10 35 (Jet 1) and 10 37 erg s _1 (Jet 2), and injection velocity 1.66 x 10 10 cm s _1 . 
The medium in which the jets propagate is an isotropic wind (as seen from the star) 
of mass loss rate 10 -6 M yr _1 and (constant) velocity 2 x 10 8 cm s _1 , typical for 
an O-type star (see Section 2 and Table 1 in PBR08). The initial jet densities are 
different, being p\ = 0.088 p a and p2 = 8.8 p a , where p a = 3 x 10 -15 gr cm -3 is 
the stellar wind density. The initial jet temperature is Tj ~ 10 10 K, Mach number 
is Mj = 16.6 and pressure P h \ = 70.8dyncm -2 and Pj 5 2 = 7.1 x 10 3 dyncm -2 . The 
star is located at a distance of 2 x 10 12 cm from the base of the jet, in a direction 
perpendicular to the jet axis. Hereafter, we will use coordinate z for the direction 
of propagation of the jet, x for the direction connecting the jet base and the star, 
and y for the direction perpendicular to both. 

We have used a finite-difference code, named Ratpenat, which solves the 3D 
equations of relativistic hydrodynamics written in conservation form. Ratpenat has 
been parallelized using a hybrid scheme with both parallel processes (MPI) and 
parallel threads (OpenMP) inside each process (see Ref. 8 for further details). The 
simulations have been performed in Mare Nostrum, at the Barcelona Supercom- 
puting Centre (BSC) using up to 128 processors. The numerical grid box expands 
transversally 20 Rj (jet radii) on each side of the jet axis (making a total of 40 Rj), 
and 10 Rj per node (4 processors). The numerical resolution of the simulation is of 
4 cells per initial jet radius. This means that the final box is 160 x 160 x 1280 cells. 



15, 2010 0:12 WSPC/INSTRUCTION FILE perucho 



On the interaction of jets with stellar winds in XRBs 3 

An extended grid is used in the transversal direction, with 80 cells, which brings 
the outer boundary 80 Rj farther from the axis. The resolution is relatively low due 
to the amount of computational time needed to perform the simulations. However, 
as the jets initially expand, the effective resolution at the distances of interest, say 
~ 10 12 from the base, is ~ 16 cells per jet radius. 

We implicitly assume that the magnetic field has no dynamical influence in the 
evolution of the jet, that the wind is continuous and homogeneous during the time 
of the simulations and that the compact object is at the same orbital position during 
this time (10 2 — 10 3 s compared to orbital times T > 10 5 s). 

2.2. Results 

Jet 1 propagates up to z ~ 1.6 x 10 12 cm after ~ 1.25 x 10 3 s. In the beginning of the 
simulation, the jet expands and generates a thick shear layer with positive velocities 
in the jet direction. The backflow surrounds and interacts strongly with this outer 
jet region, generating instabilities that grow in the shear layer. These instabilities 
are asymmetric in the plane of impact of the wind due to the different pressures in 
the cocoon on both sides of the jet (see Fig. [I]). After expansion, the central region 
of the jet becomes underpressured with respect to its surroundings and recollimates 
until the formation of a reconfinement shock. This quasi-steady shock propagates 
very slowly from z^2x 10 11 cm to z ~ 4 x 10 11 cm, as the pressure in the cocoon 
drops (see PBR08). Downstream of the reconfinement shock, the instabilities, which 
were growing in the shear layer, propagate to the whole section of the jet as the 
internal jet flow is decelerated and becomes more sensitive to perturbations. This 
process ends up in the mixing and deceleration of the jet flow at z > 10 12 cm. 

Fig. [T] shows two axial cuts of Jet 1 of rest-mass density, along the YZ (upper) 
and the XZ planes (lower) at the last snapshot. The deviation caused by the wind 
thrust can be observed in the XZ plane, which is the plane of symmetry of the wind. 
Fig. [2] shows transversal cuts of axial velocity, logarithm of rest-mass density and 
traceo at z c± 1.3 x 10 12 cm. In this figure, we can see the deformation of the bow- 
shock caused by the wind thrust. The jet has been entrained, at this axial position, 
by the wind material —fmax < 1? an d the maximum velocity in the jet fluid is still 
relatively fast (~ 1.5 x 10 10 cms _1 ), despite the irregular morphology and mixing. 
Actually, the velocity drops to v ~ 10 10 cms _1 after this point. The low velocity 
of the jet at the end of the simulation, along with its final (destabilized) structure, 
implies that the jet will be disrupted and will not propagate collimated out of the 
binary system. 

The evolution of Jet 2 is very similar to that of Jet 1 from a qualitative point of 
view, but this jet propagates up to z c± 2 x 10 12 cm in just ~ 210 s. The jet expands 
more at the base because it is initially denser (more overpressured) and the velocity 



a The tracer, / = [0, 1], indicates the composition of the fluid, with corresponding to pure wind 
material, 1 to pure jet material and any value between and 1 indicates the relative amount of 
jet material in a cell. 
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Fig. 1. Jet 1, cuts of rest-mass density (in gem -3 ) along the propagation axis. The upper panels 
show a cut perpendicular to the plane of symmetry of the wind (YZ plane), whereas the bottom 
panels show a cut in this plane (XZ plane). 
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Fig. 2. Jet 1, transversal cuts of axial velocity (in cms 1 ), rest-mass density (in gem 3 ) and 
jet mass fraction at z ~ 1.3 x 10 12 cm. 



of the jet head is faster. Therefore, the reconfinement shock is stronger and occurs 
farther downstream (see PBR08). The location of the shock changes with time from 
z^6x 10 11 cm to z ~ 10 12 cm (see Fig. [3]). The effect of the wind in the direction 
of the jet propagation is very small, as seen in Figs. [3] and 2J In the latter, the 
structure of the bow shock at z ~ 1.5 x 10 12 cm is observed to be more symmetric 
than that in Jet 1 (Fig. [2]). The devolopment of asymmetric KH instabilities in the 
shear layer, as in Jet 1, propagates to the whole jet after the reconfinment shock, 
what triggers helical motions and distortions in the jet. Fig. |4] shows that the jet 
core is unmixed (/ = 1) and that the flow velocity is still as high as that in the 
injection point. At the end of the simulation, the velocity of the bow-shock ahead of 
the jet is ~ 8.4 x 10 9 cms _1 , which is close to the initial speed (^9x 10 9 cms _1 ). 
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Fig. 3. Jet 2, cuts of rest-mass density (in gem -3 ) along the propagation axis. The upper panels 
show a cut perpendicular to the plane of symmetry of the wind, whereas the bottom panels show 
a cut in this plane. 
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Fig. 4. Jet 2, transversal cuts of axial velocity (in cms -1 ), rest-mass density (in gem -3 ) and 
jet mass fraction at z ~ 1.5 x 10 12 cm. The upper panels show a cut perpendicular to the plane of 
symmetry of the wind, whereas the bottom panels show a cut in this plane. 



3. Discussion 

Our work shows that the stellar wind from a massive star may cause the disruption 
of a jet if its power is < 10 37 ergs _1 . The cocoon generated by the jet presents 
asymmetric properties due to the impact of the wind. The difference in pressure on 
both sides of the jet triggers helical perturbation modes in the jet, which develop 
initially only in the shear layer of the jet. After the jet core goes through a strong 
reconfmement shock, the helical instability propagates to the whole jet. In PBR08, 
it was shown that such a shock will occur inside the binary region if the temperature 
of the jet is Tj < 10 14 K when the jet is surrounded by the cocoon, or Tj < 10 13 K 
when the jet is in direct contact with shocked wind. Reconfmement shocks and 
those produced in the interaction with the stellar wind (bow-shock and reverse- 
shock in the wind impact region) are candidate locations for the production of high 
energy emission. This aspect will be treated elsewhere (Bosch-Ramon, Khangulyan 
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& Perucho, in preparation). 

A candidate X-ray binary where disruption of the jet could be taking place is 
LS 5039 (see Ref. 9), but higher resolution observations are required for a proper 
probe of the jet- wind interaction region. The presented situation could take place in 
several high-mass XRBs (HMXB) in the Galaxy. The luminosity function derived 
in Ref. [TO] predicts 3 HMXBs with L x = 10 35 erg/s. Following Ref. 11, a HMXB 
with a 10 Mq black hole, could produce a jet with kinetic power between 10 35 and 
10 38 erg/s, which is in the range of the simulations performed here. Although Ref.fTOl 
does not offer any specific prediction for Lx < 10 35 erg/s, extrapolating the given 
luminosity function, we deduce that there is room for a few (~ 10) sources in our 
Galaxy in which the jets could be disrupted by the stellar wind. 

The possible influence of magnetic fields or an inhomogeneous wind remain to 
be tested. At present, we compute the evolution of Jet 2 in a decreasing density 
atmosphere given by the wind, and study the influence that this may have in the 
position of the recollimation shock and in the growth of instabilities. 
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